Persistent currents in interacting systems: role of the spin. 
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Persistent currents flowing through disordered mesoscopic rings threaded by a magnetic flux are 
investigated. Models of fermions with on-site interactions (Hubbard model) or models of spinless 
fermions with nearest neighbor interactions are considered on 2D cylinders with twisted boundary 
conditions in one direction to account for a magnetic flux. Self-consistent Hartree-Fock methods 
are used to treat the electron-electron interaction beyond first order. We show that the second 
harmonic of the current (which is relevant in the diffusive regime) is strongly suppressed by the 
interaction in the case of spinless fermions while it is significantly enhanced in the Hubbard model. 
Our data also strongly suggest that the reduction (increase) of this harmonic is related to a strong 
increase (reduction) of the spacial fluctuations of the charge density. 

PACS numbers: 72.10.-d, 71.27.+a, 72.15.Rn 



The observations of mesoscopic currents in very 
pure metalUcpiJiano-structures was done in pioneering 
experimentsu'Htl. In the first case the experiment dealt 
with the average current of a system of 10 7 disconnected 
rings in the diffusive regime while in the second a single 
ring was used. Although the existence of such persist' 
currents in small metallic rings was predicted long ag< 
the magnitude of the observed currents is still a real chal- 
lenge to theorists. There is a general belief that the inter- 
action plays a crucial role in enhancing the current. But, 
so far, the role of the interaction in disordered systems is 
still unclear. Treating interaction and disorder on equal 
footings is a difficult task. Previous work&Q has shown 
by exact diagonalizations (ED) of small clusters that, for 
strictly ID systems of spinless fermions, the effect of a re- 
pulsive interaction is to increase further the localization 
of the electrons and hence to decrease the value of t 
current. Using a Hartree-Fock approach, Kato et al. 
have obtained a qualitatively good agreement with t 
exact calculations. On the other hand, Giamarchi et al 
have pointed out that, for the ID Hubbard model, i.e. 
when spin is included, the interaction enhances the per- 
sistent current. In this case, the increase of the current is 
closely related to the decrease of the spacial fluctuations 
of the charge density or, equivalently, to the smoothing 
out of the charge density as it occurs in the ID Hubbard 
model whith repulsive interaction. This emphasizes the 
important role of the spin in ID systems. 

In higher dimensions the role of the spin is still unclear. 
First order calculations for which spin is irrelevant have 
shown thatpfche persistent currents are increased by the 
interactions^!! More recently, Ramin et al.£3 have nu- 
merically shown that the first order Hartree-Fock (HF) 
correction to the second harmonic of the persistent cut- 
rent was in agreement with the analytical treatment!!!! 
Hence, for both the spinless fermion model with nearest- 
neighbor interactions and the Hubbard model, the sec- 
ond harmonic is enhanced. However, in this treatment 



it is found that a nearest neighbour interaction tends 
to decrease the value of the typical current while a re- 
pulsive extended Hubbard interaction enhances it. In 
some previous work, Exact Diagonalisations (ED) calcu- 
lations have been compared to a Self-consistent Hartree- 
Fock (SHF) treatment of the interaction. For the small 
clusters (4x4 clusters) which could be handled we haw 
found a good agreement between the two sets of datalia. 
This direct comparison with the exact results has there- 
fore established some degree of reliability of the self- 
consistent Hartree-Fock approximation at least in the dif- 
fusive regime. In this paper, we use both HF and SHF 
treatments of the interaction between particles which en- 
ables us to treat much larger systems. The SHF treat- 
ment takes into account high order terms in the inter- 
action and, simultaneously, deals with quantum interfer- 
ence effects due to the disorder somehow exactly. It is 
important to note that this method is different from the 
usual perturbative approachli!3 where the corrections to 
the current due to the interacting term are calculated 
perturbatively. In contrast to their approach, our pro- 
cedure includes a resummation of higher order terms 
through a self-consistency relation, which turns out to 
become essential at moderate interaction. Nevertheless, 
an exact connection with some diagrammatic expension 
is a tedious problem and this issue probably deserves fur- 
ther study. Note that our study also applies in principle 
to single or multi-ring experiments. From a theoretical 
point of vue, the difference simply relies in the absence 
or presence of particle number fluctuations. In the fol- 
lowing, we shall assume that the number of electron of 
the ring is fixed. 

This paper is organized as follows: first, we compare, 
for both spinless and Hubbard models, the first order 
correction (in the interaction) and the SHF correction to 
the second harmonic of the persistent currents. We show 
that, in the case of spinless fermions, the two methods 
stay in good agreement only for rather small values of 
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the interaction parameter but a complete disagreement 
appears for moderate values. In contrast, in the Hubbard 
case, the agreement between the two approaches is rather 
good. Secondly, we present evidences that this decrease 
(increase) of the persistent currents is directly related 
to the increase (decrease) of the spacial fluctuations of 
the charge density from site to site as the interaction is 
switched on. Thirdly, the effect of the electron repulsion 
on the second harmonic is shown to increase with the 
system size. 

The Hamiltonian is defined on a L X L lattice with 
periodic boundary conditions in one direction (e.g. x 
direction) and reads: 



H — T~Lk + H 



H a 



(1) 



TLk = ^ij c i c j i s the usual kinetic part containing 
the flux dependance, ty = | exp {i-^-(i x — j x )} if i and 
j are nearest neighbor sites and otherwise. Hint is the 
interacting part and Hdes is the term due to the disorder, 
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where ni is the local density operator at site i and w\ 
are on-site energies chosen randomly in [— ^-}. When 
spin is included, n\ — n- x \ + and the kinetic term 
contains an additional sum over the spin indices. In the 
spinlcss fermion case the electron repulsion is given by 



n int 



(3) 



where Vij of strength V = |Vy| only connects nearest 
neighbor sites (screened interaction). Lastly, the Hub- 
bard interaction is defined by 



(4) 



In the SHF approximation, the interaction part of the 
spinlcss fermion model reduces to, 

n?nt = - E Sti i c i C J + E 6wi ni 

ij i 

-\j2^(^)o(^)o-\(4^o\ 2 ) (5) 

where Q Q stands for the expectation value in the ground- 
state wavefunction. The expression above contains 
two terms: an Hartree term proportional to Swi = 
Sj^i HivWo w hi° n comes out as an extra on-site dis- 
order potential and a Fock term proportional to Stij = 
T4j^CjCj) which is an extra hopping amplitude. The 

quantities (*n.j) and (cj Ci) are calculated self consis- 
tently so that the SHF Hamiltonian itself depends on 
the filling factor, the desorder strength etc... Similarly, 
in this approximation Hf nt becomes, 



H 



H 



U E(( ni T)o "U + ( Ui 0a n U - < n il)o( n iT) )- ( 6 ) 



Since we are interested in the paramagnetic phase we 
shall assume (^ix) = ( n it)o- ^ n this case the spin f and 
I are decoupled so that one can write Hf nt = J2a- ^int 
where, 



(7) 



The current is defined as the derivative of the total 
energy versus flux, 
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where E(Q) is the total energy. The current is a periodic 
function of $ of period 1 ($ is measured in units of $o) 
and thus can be expanded as a Fourier series, 



I hn sin(27rn$) 



(9) 
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where Ihn are the harmonics of the current. It is now well 
established that the ensemble average (over filling or. 
order) suppresses the first harmonic of the current 
This fact was indeed observed in multi-ring experiments 
where the current was found to be ^ periodic. Therefore, 
in the following we shall focus on the second harmonic 
Ih2- The solution of the previous set of non- linear self- 
consistent equations are obtained numerically on small 
clusters by an iterative procedure for arbitrary values of 
$ and arbitrary disorder realisations. The various physi- 
cal quantities are then averaged over disorder. The disor- 
der average denoted by Q.. in the following corresponds 
typically to an average over at least 1000 configurations 
of the disorder. 

As noted previously, the self-consistency relations can 
take care of high order effects in the interaction. It is 
therefore necessary, as a preliminary study, to investi- 
gate the role of the self-consistency by comparing the 
SHF results to the simple first order calculations (refered 
to hereafter as HF). In Fig. 0(a) (Ih2 ) d - for spinless 
particles is plotted as a function of V. The calculations 
have been done at half-filling on a 8 x 8 cylinder. The 
HF contribution is of course linear in V. As expected, we 
observe for small values of the interaction a perfect agree- 
ment between the SHF and the HF calculations. Indeed, 
the slope at V=0 of the SHF results is actually given by 
the HF data. However, for increasing V the SHF calcu- 
lation shows a strong reduction of the current whilst the 
HF predicts an increase. We also observe that the region 
of agreement between HF and SHF is reduced as the dis- 
order increases. In other words, this means that, as the 
strength of the interaction increases, the effect of higher 
order terms becomes more and more dominant and thus 
a first order calculation is not sufficient. We will see 
later on that this reduction is related, as in the ID case, 
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to an increase of the spacial fluctuations of the charge 
density. The influence of the filling factor and of the dis- 
order strength on (lh2)j- vs V is shown in Fig. |l|(b). 
This figure shows that the repulsive interaction is detri- 
mental to the persistent currents which are drastically 
suppressed for a wide range of parameters. Note that, 
at quarter filling, (lh,2) d - starts immediately to decrease 
with V contrary to the behavior observed at half-filling. 

In contrast to the spinless model, the Hubbard model 
exhibits a completely different behavior as seen in Fig. 
^ showing the relative increase of (lh2) dis as a function 
of the Hubbard repulsion U. Interestingly enough, both 
SHF and HF calculations predict an increase of the cur- 
rent. It is also interesting to notice that the first order 
calculation gives reliable results regarding the effect of 
the interaction. However, the first order calculations al- 
ways predict larger values of the currents. 

At this point, these results already suggest that the na- 
ture of the interaction plays a crutial role: in the spinless 
case the currents are strongly reduced by the interaction 
in contrast to a Hubbard repulsive interaction which en- 
hances the currents. Secondly, we have shown that, in 
the spinless fermions case, higher order terms become 
rapidly dominant even for relatively small values of the 
interaction strength. Hence a first order approach is not 
sufficient. 

Let us now try to developp a physical picture that 
could help to understand the role of the spin. We shall 
argue that the enhancement (reduction) of the persistent 
currents is related to the reduction (increase) of the fluc- 
tuations of the charge density m from site to site. One 
way of observing this effect on the charge density consists 
in plotting the distribution of the local density (n,i) - For 
that purpose, we shall consider here a 10 x 10 system and 
assume that the («i) (where the subscript k labels the 
various realisations of the disorder) are independant vari- 
ables. The distribution of the charge densities can then 
be defined as, 

1 N M . 

^)^ee%):-p) do) 

k— 1 i 

where L is the length of the system (L=10) and Ndis is 
the number of disorder configurations. As usual, the lo- 
cal densities («i) are calculated self-consistently. Note 
that, for a sufficiently large system, we expect the distri- 
bution to become independant of the choice of the disor- 
der configurations. 

For spinless fermions, we have plotted in Figs. ||(a,b) 
the distribution P(p) calculated on a 10 x 10 cylinder 
for two different fillings. We clearly observe that, with 
the interaction, both the shape of the distribution change 
and the width of the distribution increases. These effects 
are strongly emphasized in Fig. ||(b) in which two peaks 
appear, showing a tendancy towards the formation of a 
charge density wave as we approach commensurability 
(n=0.5). Note that such an instability is clearly unphysi- 



cal for a metal and could be easily removed by preventing 
perfect nesting of the Fermi surface (for example, by in- 
cluding hopping terms at larger distances). However, the 
effect of the interaction to broaden the distribution P{p) 
seems to be independant of the precise details of the band 
structure and is generic for any filling. In contrast, in the 
case of the Hubbard model (Fig. ^), we observe the op- 
posite effect: the distribution shrinks around the average 
value as the interaction is switched on. 

Let us now turn to more qualitative results. The width 
of the distribution (assuming independant variables for 
different configurations of the disorder) bp is given by 

, N dis r- 

fr-scg^ DM!-)"- <"» 

where n = jfi is the filling (N e is the number of elec- 
trons) . We have studied the effect of the interaction on 
bp as a function of V (spinless fermions) or U (Hubbard 
model) for different values of the disorder strength W 
and different fillings. In Fig. |(a) bp{V, W) is plotted 
at half-filling as a function of the parameter V (spin- 
less fermions). Here we clearly observe an increase of 
the width bp of the distribution for increasing V. Very 
crudely, this effect seems to depend only on the combined 
parameter In contrast, in the Hubbard case shown 
in Fig. ||(b), we observe a reduction of bp(U, W) as U in- 
creases. In Figs. ||(a,b) we have plotted the relative vari- 
ations of the width of the distribution Sp(V, W) /bp(0 7 W) 
and bp(U,W)/bp(0,W) as a function of the interaction 
parameters V and U for various densities. We observe 
that the effects described above become stronger for 
larger fillings i.e. when the interaction between the par- 
ticles becomes more effective. Note however that, in the 
case of the spinless fermion model, this behavior is un- 
related to the commensurability since densities n=0.25 
and n=0.4 have been chosen, clearly away from half- 
filling (n=0.5). In Figs. § we see that, in the spinless 
case, the width has increased by almost a factor 4 for 
V = 0.8 and N e = 40. In the Hubbard case the reduction 
of bp(U, W) is 25% larger compared to half-filling. Our 
present study then strongly suggests that a reliable ex- 
planation of the observed large persistent currents must 
somehow take into account the spin of the particles. 

We finish our discussion by a study of the influ- 
ence of the system size on our results. In Fig. |t] 
we have plotted the relative increase of the current 
(l h 2) d .jU,W)/(l h2 ) d jO : W) vs U at fixed density but 
for different system sizes. This figure clearly indicates 
that, as the size of the system increases, the effects in- 
duced by the interaction become stronger. We also ex- 
pect that this is also true when the connectivity of the 
lattice increases (i.e. going form 2D to real 3D rings). 
Although a systematic accurate study as a function of 
the sample size is not feasible our data for the Hubbard 
model for, let say U=0.4, are not inconsistent with the 
magnitude of the currents observed in the experiments. 
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In conclusion, we have shown that, in the spinlcss 
fcrmion model, the effect of a moderate interaction leads 
to a drastic reduction of the magnitude of the current, 
contrary to what is expected from first order calcula- 
tions. For the first time, we have established that, in 
2D, taken into account the spin degrees of freedom is 
crucial to explain the enhancement of the current due 
to the interaction, similarly to the ID case. Although 
the calculations presented in this paper deal with rather 
small 2D clusters, we expect that, for larger systems and 
higher connectivity (e.g. in 3D) and in the presence of an 
Hubbard interaction, the impurity scattering will become 
even less effective to localize the electrons. Although a 
systematic study as a function of the sample size is still 
out of reach of present day computors, our data of the 
Hubbard model at intermediate U (around 0.4) arc not 
inconsistent with the magnitude of the current observed 
in the experiments. Our study suggests that the SHF ap- 
proach provides a relatively good tool to study, on equal 
footings, the effects of the interaction and of the disorder. 
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FIG. 1. (lha) as a function of V calculated within the 

\ / dis 

SHF method for spinless fermions on a 8 x 8 cylinder. An 
average over 1000 disorder configurations has been done, (a) 
Comparison between first order HF (open symbols) and SHF 
results (full symbols) at half-filling. Circles correspond to 
W = 3 and squares to W = 4. (b) SHF results at half-filling 
(full symbols) and quarter filling (open symbols) . Circles cor- 
respond to W — 3, squares to W = 4 and triangles to W = 5. 



FIG. 2. Ratio (l h2 ) dis (U, W)/(l h2 ) dis (Q, W) as a function 
of U calculated in the Hubbard model within both HF (open 
symbols) and SHF (full symbols) methods on a 8 x 8 cylin- 
der. An average over 1000 disorder realisations have been 
performed. Circles correspond to W = 3 at quarter-filling, 
squares to W = 3 at half-filling and triangles to W — 4 at 
half-filling. 
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FIG. 3. Distribution P(p) in the case of spinless fermions. 
The calculations have been done at quarter filling n = 0.25 
(a) and n = 0.4 (b) on a 10 x 10 cylinder. We have averaged 
over 30 configurations of the disorder. The values of V and 
W are given in the figure. 

FIG. 4. Distribution P(p) in the Hubbard case. The cal- 
culations have been done at half-filling (n = 1) on a 10 x 10 
cylinder. We have considered 30 configurations of the disor- 
der. The values of U and W are given in the figure. 

FIG. 5. Effect of the interaction on the width 5p(V, W) at 
half-filling as a function of the interaction parameter (U or 
V). The calculations have been done on a 10 x 10 cylinder 
and 30 configurations of the disorder have been considered. 
The values of W are given in the figure, (a) Spinless fermion 
model; (b) Hubbard model. 

FIG. 6. Ratio Sp(V, W)/Sp(0, W) as a function of the in- 
teraction parameter (U or V) for different filling factors (the 
filling is indicated on the plot) and values of W. The calcula- 
tions have been done on a 10 x 10 cylinder with 30 configu- 
rations of the disorder, (a) Spinless fermion model (n = 0.25 
and n = 0.5); (b) Hubbard model (n = 0.5 and n = 1). 
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FIG. 7. Ratio (l h2 )(U, W) dis /(l h2 )(0, W) dis as a function 
of U (Hubbard case) calculated within the SHF method. The 
calculations have been done on 6 x 6 (circles) ,8x8 (squares) 
and 10 x 10 (triangles) cylinders at half- filling (n = 1) and for 
W = 3. 1000 realisations of the disorder have been used. 
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